2019, Instituto Mexicano de Tecnología del Agua 
, Open Access bajo la licencia CC BY-NC-SA 4.0 
Tecnología y 


CienciaszAgua (https://creativecommons.org/licenses/by-nc-sa/4.0/) 


DOI: 10.24850/j-tyca-2019-03-11 
Notas 


Estimación de la disponibilidad hídrica en la cuenca del 
lago de Zirahuén 


Modeling water availability in Zirahuén Lake Basin 


Mauricio José Ríos? 

Víctor Prado-Hernández* 
Alejandro Romero-Bautista*? 
Delfino Reyes-López* 
Fermín Pascual-Ramírez? 


lUniversidad de Ciencias y Artes de Chiapas, Tuxtla Gutiérrez, Chiapas, 
México, mauricio.joseQunicach.mx 


“Universidad Autónoma Chapingo, km. 38.5 Carretera México-Texcoco, 
Chapingo, Texcoco, Estado de México, vpradohOchapingo.mx, ORCID 
http://orcid.org/0000-0001-6045-1661 


Instituto Tecnológico del Valle de Morelia, carretera Morelia-Salamanca 
km. 6.5, Morelia, Michoacán, México, aromeroQitvallemorelia.edu.mx, 
ORCID http://orcid.org/0000-0001-9490-9650 


“Benemérita Universidad Autónoma de Puebla, Puebla de los Ángeles, 
Puebla, México, delfino_reyes2001UOyahoo.com.mx, ORCID 
http://orcid.org/0000-0002-8360-0567 


"Consejo Nacional de Ciencia y Tecnología-Universidad Nacional 
Autónoma de México, Ciudad de México, México, 
fpascualrOiies.unam.mx, ORCID http://orcid.org/0000-0002-1005-9597 


Autor para correspondencia: Fermín Pascual Ramírez, 
fpascualrOiies.unam.mx 


Resumen 


Tecnología y ciencias del agua, 10(3), 275-288. DO!: 10.24850/j-tyca-2019-03-11 275 


2019, Instituto Mexicano de Tecnología del Agua 
7 Open Access bajo la licencia CC BY-NC-SA 4.0 
Tecnología y 


CienciaszAgua (https://creativecommons.org/licenses/by-nc-sa/4.0/) 


En la gestión del agua es fundamental conocer la disponibilidad hídrica 
al interior de las unidades de planeación y manejo hidrológico. En la 
cuenca del lago de Zirahuén, para mejorar el conocimiento de las 
variables del ciclo hidrológico, se empleó el modelo de simulación 
espaciotemporal en cuencas Soil-Water-Balance (SWB), apoyado en las 
características hidrodinámicas del suelo y cobertura, en conjunto con 
información climatológica diaria. Los resultados mostraron la influencia 
del uso del suelo en los escurrimientos y la recarga de acuíferos, que se 
reflejaron en la disponibilidad hídrica, la cual mostró retrocesos 
importantes con énfasis en el periodo 2009-2015. Al analizar los 
resultados con respecto a lo reportado en 2016, se encontró que el 
Indice de Aguas Subterráneas (IDAS) ha sido subestimado en 14.64%, 
en tanto que la Disponibilidad Relativa (DR) se ha subestimado por el 
150% con respecto a los valores encontrados en este estudio. 
Finalmente, para mejorar la confiabilidad y la gestión del agua, se 
recomienda implementar estrategias de monitoreo de la cantidad, 
calidad y temporalidad del agua para verificar los resultados obtenidos. 


Palabras clave: gestión del agua, escurrimiento, agua subterránea, 
balance hídrico, uso del agua. 


Abstract 


In water management, it is essential to know the water availability 
within the hydrological planning and management units. To improve the 
knowledge of hydrological cycle variables in the Zirahuén Lake Basin we 
used the Soil-Water-Balance (SWB) spatiotemporal simulation model 
supported by the hydrodynamic characteristics of the soil and land 
cover, together with daily climatological information. Our results showed 
how land use influences runoff and recharge of aquifers reflecting in 
water availability estimates, which showed an important setback with an 
emphasis during the 2009-2015 period. When comparing our results 
with the official records reported in 2016, we found that the latter 
figures underestimated the groundwater index (GI) by 14.64% and the 
relative water availability index (RWAI) by 150% relative to the results 
of our study. We concluded that, to improve the reliability and 
management of water in the basin, it is recommend verifying our results 
by implementing strategies to monitor the quantity, quality, and 
temporality of water. 
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Introducción 


El agua y su disponibilidad constituyen la base del desarrollo de las 
sociedades, aunque su distribución adolece de relación con el 
crecimiento y expansión de la población, así como del desarrollo 
económico (Shiklomanov, 1998). La disponibilidad de agua renovable en 
México se presenta en proporción ¡inversa con el desarrollo 
socioeconómico (Conagua, 2016a), generando zonas de mucha 
demanda sin información confiable sobre el funcionamiento del ciclo del 
agua. 


El Balance Hídrico (BH) es el ejercicio de evaluación de los factores que 
intervienen en el ciclo del agua y existen métodos de cálculo 
desarrollados para diferentes escalas de observación (Sokolov 8 
Chapman, 1974; Xu 8 Singh, 1998). El modelo Soil-Water Balance 
(SWB) (Dripps € Bradbury, 2007; Westenbroek, Kelson, Dripps, Hunt, 8 
Bradbury, 2010), calcula los componentes del balance hídrico a nivel 
diario (Thornthwaite € Mather, 1957), aprovechando las ventajas de los 
Sistemas de Información Geográfica y calcula el balance hídrico espacial 
a nivel de la cuenca, requiere cinco variables: 1) precipitación y 
temperatura; 2) clasificación de uso de suelo y vegetación; 3) grupo 
hidrológico de suelo; 4) dirección de flujo, y 5) capacidad de retención 
de agua en el suelo. Realiza el cálculo de los escurrimientos utilizando el 
método de Curva Numérica del NRCS (NRCS, 1986). Se ha empleado 
con éxito en la estimación regional de la recarga de aculferos 
subterráneos en la cuenca del lago Michigan y las cuencas de Wisconsin 
(Tillman, Gangopadhyay, € Pruitt, 2016). 
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Los análisis conjuntos de los elementos del balance hídrico y la demanda 
de agua comprometida permiten conocer la disponibilidad hídrica 
(Touhami et a/l., 2014) que, en México, se estima utilizando la NOM- 
011-Conagua-2015 (Conagua, 2015), y sólo se pueden localizar de 2014 
a la fecha con ausencia de información a diferentes épocas. Para 
mejorar la gestión del agua, ante la falta de información en el tiempo, 
los modelos de simulación permiten mejorar el conocimiento del balance 
hídrico (y disponibilidad) en cuencas; por ello, el objetivo de este 
trabajo fue emplear SWB en la estimación de la disponibilidad hídrica en 
la cuenca del lago de Zirahuén a diferentes épocas y analizar los índices 
de disponibilidad con respecto a los registros oficiales reportados en 
2016. 


Materiales y métodos 


El estudio se realizó utilizando la información geoespacial de la cuenca 
del lago de Zirahuén (19% 21' 09" y 19% 29' 50" N y, 101% 29' 11" y 
101% 46' 21" O), en el estado de Michoacán de Ocampo (Figura 1); 
pertenece a la Provincia Fisiográfica Eje Neovolcánico y la subprovincia 
Neovolcánica Tarasca (INEGI, 2001). 
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Figura 1. Ubicación de la cuenca del lago de Zirahuén. 


Se tomó información de la carta topográfica a escala 1:50 000 (INEGI, 
2014) para generar un Modelo Digital 
resolución espacial de 100 m. Las coberturas de uso de suelo y 
vegetación (Tabla 1) fueron convertidas a matrices de datos de formato 
ráster ASCII a la misma resolución que el MDE, en el proceso de 
acondicionamiento de la información. Se utilizó una base de datos con 
información diaria de temperatura y precipitación, proporcionada por la 
Comisión Nacional del Agua (Conagua) convertida a unidades inglesas. 


Tabla 1. Síntesis de información utilizada. 


de Elevaciones (MDE) con 


Tipo Formato Versión Escala 
Topográfico Vectorial Serie III 1:50 000 
Edafología Vectorial Serie II 1:250 000 
Límite de cuenca Vectorial Edición 2 1:50 000 
Uso de suelo Vectorial Serie 1, II, 11, | 1:250 000 
vegetación IV, V 
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Datos climatológicos Tabular 2016 N/A 


La evaluación de la variación espaciotemporal del balance hídrico se 
realizó para los años 1972, 1980, 1992, 2000, 2003, 2009 y 2015, con 
base en la fecha de publicación de cada serie de uso de suelo y 
vegetación, y la disponibilidad de registros climatológicos diarios. La 
evapotranspiración fue estimada utilizando datos de temperatura 
máxima y mínima (Hargreaves 8 Samani, 1985), considerando la 
disponibilidad de información meteorológica en los registros diarios. 


Para calcular la disponibilidad hídrica (superficial y subterránea), se 
consultó la base de datos del Registro Público de Derechos del Agua 
(Repda) (Conagua, 216c), donde sólo se encontraron registros de 
volúmenes de descarga comprometidos del año 2000 en adelante, por lo 
que se hicieron los cálculos en periodo análogo. Para la descarga natural 
comprometida de aguas subterráneas se consideró el mismo valor (0.3 
Mm*) que en el polígono del acuífero 'Nueva Italia”, pues el área de la 
cuenca se encuentra comprendida en dicho polígono. En cuanto al 
volumen de retorno del agua superficial, se asumió un criterio 
conservador al considerar constante el valor más bajo encontrado 
(2015) para el volumen de retorno (1.45 Mm?). 


El Índice de Aguas Subterráneas (IDAS) se calculó con base en los 
resultados obtenidos de la base de datos del Repda; en el Indice de 
Disponibilidad Relativa (DR) se utilizaron los datos reportados (Conagua, 
2016b) en la cuenca, sustituyendo el valor del escurrimiento natural con 
la estimación encontrada en este trabajo. 


En general, para el cálculo de la disponibilidad hídrica y los índices de 
disponibilidad (DR, IDAS) se procedió como lo marca la NOM-011- 
Conagua-2015. Los índices de disponibilidad (DR, IDAS) fueron 
calculados sólo para el año 2015 y comparados con los resultados 
oficiales reportado a principios de 2016 (Conagua, 2016b). 


Las coberturas de uso de suelo fueron codificadas con base en la tabla 
de valores de referencia del modelo SWB (Westenbroek et a/., 2010), 
para convertirlas a formato ráster y generar mapas (Figura 2) sobre la 
dinámica del cambio de uso de suelo en los periodos de análisis. 
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SERIE |: 1978-1991 


SERIE IV: 2007-2010 


SERIE III: 2002-2003 


Proyeccion: UTM zona 14 N 
Datum: ITRF 2008 época 2010 
unidades: m 


Bosque oyamel FOREST-Coniferous 

Bosque de pino FOREST-Coniferous 

Bosque de pino-encino FOREST Mixed Deciduous/Coniferous 

Bosque de encino-pino / FOREST-Mixed Deciduous/Coniferous Simbología 

Bosque de pino-encino É FOREST-Mixed Deciduous/Coniferous 
El FOREST-Mixed Deciduous/Coniferous 


Bosque encino-pino FOREST-Broad-Leaved Deciduous Red híérica yeeGOE 
Bosque mesáfilo de montaXa É FOREST-Mixed Deciduous/Coniferous ata Ñ , [a 190 
if Y ma. 1 Cuenca de Lago de Zirahuén 101 a 
Bosque mesófilo de montaña [no =; 
Vegetación tular ] [E Cuerpo de agua ME +50 ] 240 
|GRASSLAND 1 +6 uN +0 
E-TA,RP,TA,TAP, TP! [AGRICULTURE MA ”s 


Elaboración propia 


AH, ZU URBAN/DEVELOPED-High Intensity (commercial and business) 


Figura 2. Cobertura de uso de suelo y vegetación codificados para 
SWB. 


En la capa de unidades edafológicas y puntos de perfiles de suelos Serie 
II a escala 1:250 000 (INEGI, 2004), se identificaron cuatro unidades 
edáficas: andosoles, leptosoles, luvisoles y phnaeozem. Con el cruce de la 
información de unidades edáficas y perfiles de suelo, se definieron tres 
clases texturales, una por grupo hidrológico (NRCS, 1986). Los suelos 
andosoles, de textura franco-arenosa, grupo hidrológico A con código 1; 
los leptosoles, textura franco-arcillosa-arenosa, grupo C y código 3; 
finalmente, luvisol y phaeozem, textura franco-arcillosa, grupo D, código 
4, 


Se utilizaron los valores de los componentes de la textura del suelo para 
asignar valores con base en las tablas de búsqueda por tipo de suelo 
(Thornthwaite €: Mather, 1957). A los andosoles se les asignó el valor de 
1.6; los leptosoles, 2.7, y luvisol y phaeozem, 3. 
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Resultados 


Las láminas unitarias de recarga (Figura 3) muestran que los mayores 
valores se presentaron en los años 1992 y 2000, en tanto que los años 
2009 y 2015, presentan los valores más bajos. La porción central norte 
con cobertura de bosques tiene bajos valores de recarga debido a las 
altas pendientes que registran por la presencia de las mayores 
elevaciones. 


19927 


19924 


19927" 


19924 


-101930" — -101%45' -101%2' -101%39 -101%36' -101%33' -101%30' 


N 
Simbología IN 
5 y > 
5 amu] Cuenca del Lago de Zirahuén Recarga media anual (mm) 
25 [a] 259 Proyeccion: UTM zona 14 N 
Ea 100 400 Datum: ITRF 2008 época 2010 


19924" 


unidades: m 


E «so [EJ] 600 
E >00 MI >500 


Elaboración propia 0.2456 8km 


Figura 3. Distribución espacial y temporal de la recarga. 


Tecnología y ciencias del agua, 10(3), 275-288. DO!: 10.24850/j-tyca-2019-03-11 282 


2019, Instituto Mexicano de Tecnología del Agua 
Open Access bajo la licencia CC BY-NC-SA 4.0 


¿Enologlay (https://creativecommons.org/licenses/by-nc-sa/4.0/) 


CienciaszAgua 


Los valores de precipitación de 1972 fueron similares a 2015, pero los 
valores de recarga son del orden de tres veces menor en 2015 que en 
1972 (Tabla 2), mostrando los efectos del cambio de uso de suelo. 


Tabla 2. Componentes del balance hídrico en lámina unitaria mm/año. 


Año Precipitación | Recarga | Intercepción Escurrimiento ETP 

potencial 
1972 999.00 178.82 26.92 484,12 1238.25 
1980 1219.60 156.46 27.18 596.90 1771.65 
1992 1922.10 415.80 32.51 935.23 1 560.07 
2000 1 406.00 206.76 31.24 681.74 1 695.70 
2003 1071.80 107.70 44.45 519.68 1555.50 
2009 958.89 53.09 43.18 465.58 1 613.15 
2015 999.50 45.21 41.40 487,43 1 636.27 


La disponibilidad hídrica en la cuenca (Tabla 3) es el resultado del 
análisis minucioso de la base de datos del REPDA. En la demanda de 
agua subterránea se observa incremento entre el periodo de 2009 a 
2015 por el 284%, motivado por el aumento en los volúmenes 
requeridos por la agricultura, que en el periodo pasó de 0.1033 a 0.3975 
Mm?, En cuanto a la demanda de agua superficial, el abastecimiento 
público registró incrementos de 40% entre 2009 y 2015. Estos 
resultados, muestran que se está intensificando la presión sobre el agua 
al interior de la cuenca. 


Tabla 3. Disponibilidad hídrica en millones de m*/año. 


Año R Di DHi E D2 DH2 

2000 56.52 0.10 56.12 186.36 0.75 187.06 
2003 29.44 0.19 28.95 142.06 2.19 141.32 
2009 14.51 0.18 14.03 127.27 2.44 126.29 
2015 12.36 0.75 11.31 133.25 3.18 131.52 


Nota: R, Recarga; D1, demanda/extracción subterránea; DH1, 
disponibilidad hídrica subterránea; E, escurrimiento; D2, 
demanda/extracción superficial; DH2, disponibilidad hídrica superficial. 
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Discusión 


Las diferencias en las láminas de recarga se explican por el cambio de 
uso de suelo (Figura 2) en el área de bosques, pastizales y agricultura, 
que afectaron la tasa de recarga. Se ha encontrado que los cambios de 
uso forestal hacia la agricultura afectan la conservación de suelo y agua, 
así como en la disminución de las tasas de infiltración (Sun et al., 
2018). En cuencas degradadas, como medida correctiva, la rotación de 
cultivos reduce la erosión hasta en 69% (Abdelwahab, Bingner, Milillo, € 
Gentile, 2014); por otro lado, la conversión de cobertura forestal del 
suelo a tierras de cultivo (agrícola o de pastoreo) produce reducciones 
en el carbono orgánico en suelos minerales y carbono en las raíces finas 
(Mendoza-Vega, Karltun, € Olsson, 2003). El cambio de uso de suelo en 
escenarios de cambio climático disminuye en 19% la descarga anual de 
escurrimientos y el flujo máximo en el verano hasta de 45% (Pijl, 
Brauer, Sofia, Teuling, € Tarolli, 2018). 


El acuífero “Nueva Italia” tiene un IDAS de 0.82 y se clasifica en la zona 
de disponibilidad 4, con suficiente volumen para satisfacer las demandas 
de agua comprometidas en 2016 (Conagua, 2016b). En los cálculos 
utilizando los resultados de la simulación (2015) se encontró IDAS de 
0.94, 14.16% más que lo calculado en las cifras oficiales para todo el 
acuífero. Analizando los volúmenes concesionados en el periodo 2000- 
2015 (Tabla 3) hay tendencia lineal con pendiente negativa que indica 
disminución en la cantidad de agua disponible del acuífero subterráneo 
en el tiempo. 


Respecto al agua superficial, en los resultados de la simulación se 
observa tendencia decreciente, con un mínimo en 2009; a pesar de la 
aparente recuperación en 2015, los valores están por debajo de lo que 
aconteció en 2000 y 2003. Las cifras oficiales indican DR en 0.92, en 
zona de disponibilidad 1; sustituyendo los volúmenes de escurrimiento 
de agua obtenidos con SWB (Tabla 3), se obtuvo la DR en 2.32; esto es, 
150% más que en las cifras oficiales, para ubicarse en la zona de 
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disponibilidad 2, por lo que tiene volumen suficiente para satisfacer las 
necesidades actuales y un margen para las necesidades futuras. 


Los valores de escurrimientos aquí calculados (133.246 Mm* en 2015), 
son muy distantes de los considerados por Conagua (52.16 Mm”) en la 
determinación de la disponibilidad, y se ha subestimado la capacidad de 
la cuenca, por lo que la disponibilidad de agua es de 78.755 Mm”, Estas 
diferencias se deben principalmente al método utilizado para la 
estimación del escurrimiento; Conagua utiliza el método racional con 
resultados a nivel anual; SWB emplea el método de curva numérica en 
una solución geoespacial, considerando las particularidades de la 
topografía, así como el tipo y la cobertura del suelo, con lo que obtiene 
resultados del balance hídrico diario. Las diferencias en la temporalidad 
del cálculo permiten mayor sensibilidad en los resultados. En ambos 
procesos se obtiene información sobre la cantidad, dejando de lado las 
consideraciones en cuanto a la calidad del agua. 


A pesar de contar con una modelación más fina, los alcances de este 
trabajo no se han verificado con aforos, pues hasta la fecha existe un 
vacío de información al respecto; de igual manera, la información sobre 
calidad del agua y distribución temporal es deficiente. SWB ha probado 
su utilidad en la simulación del ciclo hidrológico con resultados 
importantes en Estados Unidos de Norteamérica (Tillman et a/., 2016; 
Zhang, Felzer, € Troy, 2016), con variaciones de +33 mm al año en la 
estimación de la recarga; se ha encontrado alta confianza en los 
resultados de la modelación, a escala regional, cuando se utilizan datos 
a resolución mensual y anual (Han et a/., 2018). Los estudios sobre la 
aplicación del modelo SWB dan cuenta de su utilidad como herramienta 
para apoyar los procesos de gestión del agua en cuencas. 


Conclusiones 


El uso del modelo SWB permitió realizar estimaciones del balance hídrico 
con buena aproximación en los flujos subterráneos y distantes en el 
flujo superficial del agua. Los análisis con respecto a la disponibilidad 
hídrica permitieron ajustar los valores de IDAS en 14.92% por encima 
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de las cifras oficiales y la DR en 150% más sobre lo oficial, como 
resultado de considerar las características biofísicas de la cuenca y la 
precipitación a resolución temporal diaria. Aunque la literatura reporta la 
utilidad de SWB en la estimación de la recarga, los resultados mostraron 
la urgente necesidad de implementar en la cuenca del lago de Zirahuén 
instrumentos y/o estrategias para el monitoreo de cantidad, calidad y 
temporalidad del agua para facilitar el proceso de gestión. 
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